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A  COMPARATIVE  STUDY  OF  LINEAR  ARRAY  SYNTHESIS 
TECHNIQUE  USING  A  PERSONAL  COMPUTER 


INTRODUCTION 


Procedures  for  synthesizing  the  radiation  patterns  of  linear  arrays 
based  on  the  specification  of  its  sidelobe  structure  are  well  establish¬ 
ed.  One  of  these  is  the  technique  originally  proposed  by  Dolph*  and 
it  provides  for  a  uniform  sidelobe  level.  Another  technique,  although 
developed  for  line  sources,  is  due  to  Taylor?  and  can  be  adopted  for 
linear  arrays.  The  latter  is  popular  due  to  its  synthesized  aperture 
distributions  which  are  more  readily  realizable.  A  discrete  version  of 
the  Taylor  synthesis  procedure  is  discussed  by  Shelton^. 

Both  Dol ph-Chebyshev  and  Taylor  synthesis  techniques,  fundamental¬ 
ly,  rely  on  manipulation  of  the  zeros  of  the  linear  array  pattern 
function.  The  aperture  distribution  for  the  desired  pattern  function 
usually  requires  lengthy  computations.  In  case  of  Dol  ph-Chebyshev 
synthesis,  this  problem  has  been  addressed  by  several  authors4-^ 
over  the  past  few  decades.  The  Taylor  synthesis,  in  effect,  uses  a 
discrete  Fourier  Transform  technique  (called  WoodwardlO  synthesis) 
to  obtain  the  aperture  distribution.  These  procedures  do  not  have 
much  in  common;  as  a  matter  of  fact,  in  case  of  an  endfire  Chebyshev 
array,  expression  for  the  element  excitations  are  quite  different  from 
that  for  a  broadside  Chebyshev  array.  However,  the  knowledge  of  the 
pattern  null  locations  in  the  above  synthesis  procedures  can  be  used 
to  develop  a  simple  expression  that  is  suitable  for  all  cases.  The 
expression  is  readily  developed  based  on  the  convolution  synthesis 
procedure  discussed  by  Laxpati H  for  planar  arrays. 

With  several  alternate  expressions  being  available  for  the  aperture 
distribution  of  Dol ph-Chebyshev  and  Taylor  syntheses,  it  is  desirable 
to  undertake  a  study  to  make  some  recommendations  as  to  the  suitability 
of  these  expressions  in  numerical  computation.  Due  to  the  increasing 
use  of  personal  computers  by  antenna  engineers,  it  is  felt  that  an  in¬ 
vestigation  of  this  nature  should  be  confined  to  the  computation  using 
such  small  computers.  Thus,  in  this  paper,  we  present  the  results  of  a 
comparative  study  of  various  linear  array  synthesis  techniques.  In  the 
next  section,  after  a  brief  discussion  of  the  three  basic  techniques 
for  evaluation  of  Chebyshev  coefficients,  we  discuss  the  accuracy  and 
computation  times  associated  with  these  techniques.  The  following 
section  presents  the  results  of  the  study  involving  two  different 
techniques  (one  due  to  Shelton^  and  the  other  using  the  convolution 
procedure)  for  Taylor  synthesis.  In  the  last  section,  some  general 
observations  about  the  investigation  and  on  the  results  are  offered. 

DQLPH-CHEBYSHEV  SYNTHESIS 


Following  Dolph's  paper  on  Chebyshev  synthesis,  Barbiere4, 
Van  Her  Maas$,  SalzevG,  and  Brown? .8  reported  on  alternative  means  of 
evaluating  aperture  distribution  for  Chebyshev  arrays.  Although  they 
are  not  the  same,  the  expressions  by  Barbiere,  Salzev  and  Brown  are 
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similar  in  that  they  express  the  current  in  an  element  in  terms  of  a 
finite  series  of  terms  involving  ratios  of  factorial  functions  and  with 
alternating  sign.  The  expression  by  El  1 i ott 12  is  representati ve  of 
this  group  and  is  the  one  used  in  this  work  and  is  reproduced  below. 
We  shall  call  this  the  classical  expression.  Also,  although  our  results 
are  valid  for  odd  or  even  number  of  elements,  for  simplicity,  we  will 
present  examples  of  odd  number  of  elements.  Thus,  all  linear  arrays 
discussed  in  the  following  have  (2N+1)  elements;  the  element  numbering 
scheme  is  shown  in  figure  1,  where  the  elements  are  assumed  to  have  a 
symmetric  excitation  leading  to  the  broadside  radiation. 

Classical  Technique: 


N 

v  fnN-P-iL 

L  (-1)  N+p 
p=n 


r(N+p+l)  _ _  (uo) 

r(N-p+l )r(p+n+l )f(p-n+l) 


2p 


(1) 


n=0 ,1 ,2,. . . ,N. 

where  T2n(u0)=R  and  SLL  =  20  log  R.  Here,  SLL  is  the  desired  sidelobe 
level  in  dB,  T2n(x)  is  the  Chebyshev  polynomial  of  degree  2N  and  r(x) 
is  the  Gamma  function. 

In  contrast,  the  expression  given  by  Van  Der  Maas  involves  terms 
of  the  same  sign  inside  the  summation.  Rresler9  reformulated  the 
expression  into  a  recursive  form  using  nested  products.  This,  we  feel 
is  a  distinctly  different  form  of  representation  of  the  coefficients. 
Thus,  we  use  this  representation  (called  Nested  Product  Technique)  in 
our  comparison.  This  expression  (in  our  notation)  is  shown  below. 

Nested  Product  Technique: 

IN-n  =  2N  a  NP(n,fm,a).  n=0,l,2,  (2) 


n 


n 


where  NP(n,fm,a)  =  \  an'm  it  fj;  s  1. 

m=l  j=m 


.  _  m(2N-2n+m) 

and  fm  =  (n-m) (n+l-m)  5 


also  a  =  1  - 


The  third  technique  is  based  on  the  convolution  of  three  element 
canonical  arraysll.  These  canonical  arrays  have  outer  element  excita¬ 
tions  of  unity,  whereas  the  center  element  excitation  Cj;  for  j=l,2, 
...N  is  chosen  such  that  the  jth  canonical  array  has  a  pattern  null 
at  the  location  of  the  jth  symmetric  zero  pair  of  the  Chebyshev 
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polynomial.  These  arrays  are  then  convolved  to  generate  the  large 
array.  Convolution  Technique: 

( 2  J  - 1 ) 

w0j  =  cos  “7*1  IT  ;  j=l  ,2 , . . .  ,N . 

Where  w0j  are  the  zeros  of  the  Chebyshev  polynomial  T£n(w). 

Cj  =  -2cos  uQj ; 

u0J-  =  2arccos  (w0j/u0). 

And  the  aperture  distribution 


N 

IU)  -  1  I„«(x-nd)  =  *fN  (3) 
n=-N 

where  fj  =  o(x-d)+Cjo(x)  +  5(x+d). 

Using  these  three  expressions  (equations  (1),  (2)  and  (3)),  comput¬ 
er  programs  NESTED,  CHEB  and  CONCHEB,  respectively,  were  written  to 
implement  the  Chebyshev  synthesis.  Different  versions  of  the  program 
suitable  for  implementation  on  different  machines  were  written.  These 
were  two  personal  computers  used  in  the  numerical  phase;  one  is  an  8- 
bit  Radio  Shack  TRS-80  Model  II  which  has  available  an  interpretive 
RASIC  language.  The  other  computer  is  a  16-bit  NEC  Advanced  Personal 
Computer  with  BASIC  and  FORTRAN  IV  compilers.  Also,  in  order  to 
ascertain  the  numerical  accuracy,  some  of  the  programs  were  run  on  a  32- 
bit  mainframe  computer  (Texas  Instrument's  Advanced  Scientific  Computer 
at  the  Naval  Research  Laboratory)  using  double  precision  (REAL*8) 
arithmetic. 

The  computation  was  carried  out  for  several  different  array  sizes 
ranging  from  15  to  99  elements;  although,  in  principle,  there  is  no 
limit  to  the  size  of  arrays  that  may  be  synthesized.  Furthermore,  all 
designs  specified  a  sidelobe  level  of  30  dR. 

Figure  2  shows  the  run  time,  under  FORTRAN,  for  the  three  afore¬ 
mentioned  Chebyshev  synthesis  programs  versus  number  of  elements.  The 
CHEB  program  was  the  slowest;  but  more  importantly,  the  program  failed 
to  converge  to  the  correct  element  excitations  beyond  30  elements. 
Over  21  elements  the  accuracy  of  the  excitation  was  only  to  2  digits. 
When  the  program  was  run  using  double  precision  arithmetic  it  still 
failed  to  converge  above  31  elements.  This  indicates  that  the  classical 
technique  inherently  has  a  limitation  as  to  the  largest  size  of  array 
that  may  be  synthesized. 

The  convolution  synthesis  program,  CONCHEB,  although  much  faster 
than  CHEB,  certainly  cannot  complete  with  NESTED  program  in  speed. 
Also,  beyond  61  elements,  the  CONCHEB  program  failed  to  converge.  The 
current  version  of  the  program  convolves  three-element  arrays  using 
zeros  of  the  Chebyshev  polynomial  in  an  alternating  sequence;  i.e.,  the 
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Figure  2  -  Run  Times  for  Chebyshev  Programs 
FORTRAN  IV  on  NEC-APC 
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sequence  in  j  is  1,N,2,N-1,3,  — .  No  attempt  was  made  to  modify  this 
convolution  process  to  improve  the  accuracy;  past  experience  with  the 
convolution  process  indicates  that  some  improvement  may  be  possible. 
However,  with  reference  to  figure  2,  it  is  obvious  that  the  NESTED 
program  is  the  most  efficient  one. 

The  results  of  the  element  excitations  indicate  that  't  is  an 
extremely  stable  algorithm;  provides  a  very  good  accuracy  in  single 
precision  (six  digit  accuracy);  and  of  course,  it  is  very  fast.  This 
program,  NESTED,  was  translated  into  BASIC  and  run  on  the  TRS-80,  Model 
II  computer.  The  execution  time  ranged  from  two  seconds  for  21  element 
array  to  36  seconds  for  a  99  element  array.  Although,  the  execution 
times  in  BASIC  are  about  15  to  20  times  longer  than  that  in  FORTRAN, 
they  are  not  significantly  long  to  be  of  any  major  consequence.  The 
NESTED  program  was  also  run  using  double  precision  (16  significant 
digits)  on  the  mainframe  computer.  The  total  execution  time  for  all 
10  different  arrays  was  less  than  0.3  seconds! 

Our  experience  with  synthesis  of  various  Chebyshev  arrays  using 
these  three  different  techniques  clearly  demonstrates  that  the  most 
important  consideration  on  small  computers  is  not  the  speed  of  execution 
but  the  accuracy  of  the  final  result.  In  this  sense  as  well,  the  nested 
product  algorithm  proposed  by  Bresler^  is  the  winner. 

TAYLOR  SYNTHESIS 

Synthesis  procedure  proposed  by  Taylor2  applies  to  a  continuous 
aperture.  In  practice,  this  procedure  is  used  for  discrete  aperture 
(arrays)  by  properly  discretizing  the  continuous  distribution.  Shelton3 
presented  a  synthesis  procedure  for  discrete  aperture  distribution  for 
Taylor  type  sidelobe  structure.  He  expressed  the  pattern  function  in 
the  form  of  a  product  function  of  zeros  and  then  carried  out  the  synthe¬ 
sis  exactly  analogous  to  that  by  Taylor;  that  is,  to  use  the  Woodward 
synthesis  technique.  In  particular,  for  a  2N+1  element  array,  all  2N 
zeros  are  explicitly  specified  in  the  pattern  function.  Thus,  analogous 
to  the  Chebyshev  synthesis,  this  synthesis  is  amenable  to  the  convolu¬ 
tion  procedure.  In  view  of  this,  in  the  case  of  Taylor  synthesis,  we 
compare  the  two  techniques;  one  proposed  by  Shelton  and  the  other 
being  the  convolution  synthesis.  Before  presenting  and  discussing  the 
results  of  the  investigation,  the  pertinent  expressions  for  the  two 
syntheses  are  given  below.  Once  again,  we  will  limit  out  discussion 
to  arrays  with  odd  (2N+1)  number  of  elements. 

Discrete  Taylor  (Shelton-^)  Technique: 

2*77  ^  A2+(n-l/2)2 

uon  =  T2N+1)  j- - — 

\  A2+(n-l/2)2 

2-rrn 

=  (2N+1) 


,  n  =  l ,2 ,--- ,n-l 


(4) 


,  n=n,  —  ,N. 


f 


J. 

where  A  =  w  cosh-l(R);  n  is  equal  to  the  number  of  near-in  zeros  that 
are  moved  in  order  to  achieve  the  desired  sidelobe  ratio  R  (or  equiva¬ 
lently  the  number  of  near-in  sidelobes  that  are  required  at  the  speci¬ 
fied  level).  The  element  excitations  are 


n-1 


2mpir 


Ip  =  1+2  l  am  cos  2N+1  >  p=0,  — ,N 


m=l 


(5) 


where 


.  2itm 
am  =  E  l  2N+1  - 


N 

E(u)  =  n 
n=l 


(cosu-cosuon) 

(l-cosuon) 


For  the  case  of  the  convolution  synthesis  procedure,  once  the 
symmetric  zero  pairs  are  established,  the  excitation  of  the  center 
element  of  a  three  element  canonical  array  is  readily  determined.  The 
procedure  and  expressions  are  analogous  to  the  case  of  Chebyshev 
convolution  synthesis.  They  are 

zeros  arej+u0j-,  j=l,2,---,N 

» 

where  u0j  are  defined  through  equation  (4),  and  the  excitation 

Cj  =  -2cosu0j. 

The  synthesis  of  the  large  array  is  carried  out  using  the  convolu¬ 
tion  of  three  element  arrays,  chosen  in  the  same  alternating  zero 
sequence  as  indicated  for  the  Chebyshev  array. 

Based  on  these  two  procedures,  computer  codes  STAYL  and  CONTAYL, 
respectively,  were  developed  in  FORTRAN  using  single  precision  arithme¬ 
tic  Run  time  associated  with  these  codes  for  “  =  6  and  the  sidelobe  level 
of  30  dB  for  various  number  of  elements  from  15  to  99  were  recorded  and 
are  shown  in  figure  3. 

The  program  CONTAYL  failed  to  converge,  once  again,  for  arrays 
with  more  than  71  elements  and  provided  only  two  to  three  digit  accuracy 
between  31  and  61  elements.  These  results  are  similar  to  the  Chebyshev 
convolution  synthesis.  Even  the  run  time  data  is  very  close. 

The  computation  time  associated  with  STAYL  has  an  interesting 
behavior  with  increasing  number  of  elements;  it  is  almost  linear.  This 
is  to  be  expected,  since  the  number  of  computations  to  be  carried  out 
for  each  element  is  determined  by  "n  and  not  (2N+1).  The  corresponding 
growth  for  CONTAYL  is  exponential.  Thus,  for  small  number  of  elements 
CONTAYL  may  save  some  computation  time  but  will  suffer  in  accuracy  as 
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Figure  3  -  Run  Times  for  Taylor  Programs  - 
FORTRAN  IV  on  NEC-APC 
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the  number  of  elements  increases.  A  check  of  STAYL  program  using  double 
precision  arithmetic  on  the  mainframe  computer  indicates  that  it  has 
five  to  six  digit  accuracy  in  single  precision  on  a  small  computer. 

It  should  be  noted  that  the  STAYL  program  code  was  developed  by 
Shelton  for  the  HP-41C,  a  pocket  calculator.  On  this  calculator,  one 
has  10  significant  digit  capacity  and  thus  the  results  obtained  are 
more  accurate  than  with  a  single  precision  FORTRAN  Code.  But,  as  one 
would  expect,  the  HP-41C  is  very  slow;  it  took  approximately  5  minutes 
to  synthesize  a  31  element  array. 

STAYL  Code  was  also  run  on  NEC-APC  using  CBASIC,  a  compiler  BASIC. 
In  CBASIC,  the  computation  times  were  significantly  higher,  ranging 
from  30  seconds  for  a  15  element  array  to  217  seconds  for  a  99  element 
array.  However,  the  computation  was  carried  out  to  14  significant 
fi gures . 

Once  again,  as  with  Chebyshev  synthesis,  we  find  the  overriding 
consideration  in  Taylor  synthesis  is  not  the  computation  time,  but  the 
accuracy  of  the  results.  In  this  sense,  Shelton's  procedure  is  most 
efficient. 

CONCLUSIONS 

As  is  often  the  case  with  engineering  investigations,  the  most 
significant  results  presented  in  this  paper  are  not  what  we  were  locking 
for  when  we  began  the  project.  We  were  originally  interested  in 
evaluating  computer  run  times  for  the  various  programs.  However,  two 
points  soon  became  apparent  --  first,  most  of  the  programs  run  fast 
enough,  even  on  small  machines,  so  that  run  time  is  not  a  major  concern, 
and  second,  only  two  of  the  programs  give  adequate  precision  for  the 
range  of  array  size  that  was  investigated.  It  is  concluded  that 
Bresler's  nested  product  algorithm  gives  excellent  results  in  terms  of 
speed  and  precision,  and  also  that  Shelton's  discretized  procedure 
allows  precise  Taylor  synthesis  for  all  sizes  of  arrays.  Finally,  it 
is  noted  that  the  programs  are  very  brief;  the  FORTRAN  computer  codes 
for  all  five  programs  are  included  in  the  appendix  and  the  codes  in 
BASIC  are  also  available  from  the  authors. 
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APPENDIX 


In  this  appendix  the  FORTRAN  IV  computer  codes  for  NEC-APC  with 
supersoft  FORTRAN  compiler  are  listed.  As  noted  in  the  main  body  of 
the  report,  the  programs  are  brief;  there  are  a  number  of  "comment" 
statements  in  the  listing  and  thus  are  easy  to  follow.  No  sample  inputs 
or  outputs  are  included. 

Programs  listed  in  the  following  are  NESTED,  CHEB,  CONCHEB,  STAYL 
and  CONTAYL. 


11 


■ 


PROGRAM  NESTED 


001 

002 

002 

003 

00-4 

005 

006 

007 

008 

009 

010 

Oil 

012 

013 

014 

015 

016 

017 

018 

019 

020 

021 

022 

023 

024 

025 

026 

027 

028 

029 

030 

031 

032 

033 

034 

035 

036 

037 

038 

039 

040 


C 

C  REVISED  01/28/84 
C  CHEBYSKEV  ARRAY  ALGORITHM  USING  NESTED  PRODUCTS  FORMULATION 
C  INPUT  IS  M  =  NO.  OF  ELEMENTS;  SLL  =  SIDELOBE  LEVEL  IN  DB. 
REAL  1(100),  NP,  C(100) 

WRITE  (1,100) 

100  FORMAT  (’  ENTER  DATA:  M.SLL' ) 

READ  (1,200)  M,SLL 
200  FORMAT  (  IO.FO.O) 

N  =  M  /  2 
TEST  =  (-1 )**M 
IF(TEST.GT.O)  GO  TO  10 
N  =  (M  -  1)  /  2 
10  R  =  10.  ••  (  SLL  /  20.  ) 

ARCOSH  =  ALOG  (  R  SQRT  (  R*«2.  -  1  )  ) 

A  =  ARCOSH  /  (M-1 ) 

ALPHA  =  (  TANH  (A)  )  **  2. 

I(N+1)  =  1.0 
I(N)  =  (M-1)  •  ALPHA 
DO  30  K  =  2,  N 
NP  =  1.0 

DO  20  J  =  1 ,  K-1 

FN  =  J  *  (M-1-2*K+J) 

FD  =  (K-J)  *  ( K-f  1  -  J ) 

F  =  FN/FD 

NP  =  NP  *  ALPHA  «  F  +  1. 

20  CONTINUE 

I(N+1-K)  =  (  M-1  )  »  ALPHA  •  NP 
30  CONTINUE 

DO  40  L  =  1,  N+1 

C(N+L)  =  I(L) 

C(N+2-L)  =  I(L) 

40  CONTINUE 

WRITE  (4,  50) 

50  FORMAT  (  '  CURRENTS’) 

DO  60  L  =  1,  M 

WRITE  (4,70)  L,  C(L) 

60  CONTINUE 

70  FORMAT  (10X,  12,  10X,  F10.6) 

STOP 

END 
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001 

c 

PROGRAM  CHEB 

002 

C  REVISED  01/28/84 

003 

C  BASED  ON  A  CLASSIC  METHOD  OF  COMPUTATION  OF  I 

004 

C  EXCITATION  VOLTAGES. 

005 

C  REFERENCE  ANTENNA  THEORY  AND  DESIGN;  ELLIOTT 

006 

C  ODD 

NUMBER  OF  ELEMENTS  ONLY 

007 

REAL  CC(100),C(100) ,CRNT( 100) 

008 

WRITE  (1,100) 

009 

100 

FORMAT  (’  ENTER  DATA:  N,  SLL* ) 

010 

READ  (1,200)  N,SLL 

011 

200 

FORMAT  (IO,FO.O) 

012 

M  =  N-1 

013 

MM  =  M/2 

014 

NN  =  (N+1)/2 

015 

PI  =  3-1415927 

016 

R  =  10.  »*  (SLL/20.0) 

017 

U  =  COSH  ( RCOSH  (R)/  FLOAT  (M)  ) 

018 

DO  10  I  =  1 ,NN 

019 

II  =  I  -1 

020 

o 

. — s 

M 

II 

O 

O 

021 

DO  20  J  =  I,NN 

022 

JJ  =  J  -  1 

023 

A  =  FLOAT  (NN  +  JJ) 

024 

GA  =  GAMALN  (NN  +  JJ) 

025 

GB  =  GAMALN  (NN  -  JJ) 

026 

GE  =  GAMALN  (J  -  II) 

027 

GD  s  GAMALN  (J  +  II) 

028 

UP  =  U  **  (2*J J) 

029 

SIGN  =  (-1)  *»  (NN-J) 

030 

TL  s  EXP  (GA  -  GB  -GE  -  GD) 

031 

TN  =  UP*SIGN*(2.*NN-1 . )/(2.*A) 

032 

T  =  TL  •  TN 

033 

C(I)  =  C(I)  +  T 

034 

20 

CONTINUE 

035 

10 

CONTINUE 

036 

DO  12  J  =  1  ,NN 

037 

CC(J)  =  C(J)/  C(NN) 

038 

12 

CONTINUE 

039 

DO  13  J  =  1,NN 

040 

CRNT  (NN-1+J)  =  CC(J) 

041 

CRNT  (NN+1-J)  =  CC(J) 

043 

13 

CONTINUE 

044 

WRITE  (4,30) 

045 

30 

FORMAT  (•  CURRENTS’) 

046 

DO  14  I  =  1,N 

047 

WRITE  (4,300)  I,  CRNT(I) 

049 

14 

CONTINUE 

050 

300 

FORMAT  (10X,I2, 15X.F12.8) 

051 

STOP 

052 

END 

CHEBYSHEV 
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PROGRAM  CONCHEB 


001 

PROGRAM  CONCHEB 

002 

!«««« 

003 

c* 

LINEAR  ARRAY  SYNTHESIS  USING  CONVOLUTION  METHOD. 

004 

c* 

CHEBYSHEV  SIDELOBE  DESIGN 

005 

c* 

ODD  NUMBER  OF  ELEMENTS. 

006 

!«*«* 

007 

REAL  PSI(IOO) 

008 

REAL  C(1 00 ), AA (1 00 ),A1(1 00 ),A2 (100), A3 (100) ,CONV(  100) 

009 

DATA  A1/100*1 ,/ ,  A3/100*1 ./ ,  CONV/100*0./ ,  AA/100*0./ 

010 

C»n 

I«l»« 

Oil 

c» 

N  =  NUMBER  OF  ELEMENTS  IN  THE  ARRAY.  MUST  BE  ODD!! 

012 

c* 

SLL  =  SIDE  LOBE  LEVEL  IN  DBS. 

013 

014 

PI  =  3.1415297 

015 

WRITE  (1,100) 

016 

100 

FORMAT  ('  ENTER  DATA:  N,SLL») 

017 

READ  (1,200)  N,SLL 

018 

200 

FORMAT  (IO,FO.O) 

019 

M  =  N-1 

020 

NN  =  (N+1)/2 

021 

MM  =  (N-1)/2 

022 

MD  =  (MM/2)  +  1 

023 

20 

CALL  CHEBX(PI,M,NN,SLL,MM,PSI) 

024 

025 

c* 

CREATE  THREE  ELEMENT  ARRAYS 

026 

»«•*«« 

027 

30 

DO  40-1  =  1 ,MM 

028 

A2(I)  =  -2.  •  COS(PSI(D) 

029 

40 

CONTINUE 

030 

IHIfl 

031 

C» 

REPEATED  CONVOLUTION  OF  3-ELEMENT  ARRAYS 

032 

c*»* 

>»••»* 

033 

CONV(1)  =  A1(1) 

034 

C0NV(2)  =  A2( 1 ) 

035 

C0NV( 3)  =  A3 ( 1 ) 

036 

L  =  1 

037 

K  =  5 

038 

LX  =  0 

039 

50 

LL  =  NN-L 

040 

LX  =  LX+1 

041 

60 

L  a  LL 

042 

C(1)  =  A1(L) 

043 

C(2)  =  A2(L) 

044 

C(3)  =  A3(L) 

045 

DO  70  I  =  1,3 

046 

DO  70  J  =  I,K 

047 

JJ  =  J-I+1 

048 

AA( J)  =  AA( J)  +  CONV(JJ)*C(I) 

049 

70 

CONTINUE 

050 

DO  80  I  =  1,K 

051 

CONV(I)  =  AA(I) 

052 

AA(I)  s  0.0 

053 

80 

CONTINUE 

054 

K  s  K+2 

055 

IF  (L.EQ.MD)  GO  TO  90 

056 

LL  =  NN+1-L 

057 

LSUM  =  L+LX 

059 

IF  (LSUM.EQ.NN)  GO  TO  60 

060 

GO  TO  50 

061 

90 

CONTINUE 

062 

WRITE  (4,600) 

063 

600 

FORMAT  (  »  CURRENTS'/) 

064 

DO  120  I  =  1,N 

065 

WRITE  (4,700)  I,CONV(I) 

066 

120 

CONTINUE 

067 

700 

FORMAT  ( 1 0X, 12, 1 0X, FI 0.6) 

068 

STOP 

069 

END 

070 

FUNCTION  COSH(R) 

071 

Y  =  EXP(R) 

072 

COSH  =  (Y  +  ( 1 .0/Y) )/2. 

073 

RETURN 

074 

END 

075 

c 

r««m 

076 

c» 

INVERSE  HYPERBOLIC  COSINE  FUNCTION 

077 

IHHi 

078 

FUNCTION  RCOSH(R) 

079 

RCOSH  =  ALOG( R  +  SQRT(R*R  -  1.0)) 

080 

RETURN 

081 

END 

082 

c*** 

083 

c* 

CHEBYSHEV  ZEROS 

084 

085 

SUBROUTINE  CHEBX ( P I , M , NN , SLL , MM , PSI ) 

086 

REAL  X(50) ,PSI( 100) 

087 

R  =  10.0  **(SLL/20. ) 

088 

B  =  COSH( RCOSH( R)/M) 

089 

DO  10  I  =  1 ,NN 

090 

J  =  1-1 

091 

X(I)  =  C0S(PI*(2.*J  ♦  1 . )/(2.*M) ) 

092 

10 

CONTINUE 

093 

DO  20  J  r  1 ,MM 

094 

II  =  NN-1+J 

095 

JJ  =  NN-J 

096 

Y  =  X(J)  /  B 

097 

PSI(II)  =  2.*ATAN(SQRT( 1-Y*Y)/Y) 

098 

PSI(JJ)  =  PSI(II) 

099 

20 

CONTINUE 

100 

RETURN 

101 

END 
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001 

c 

PROGRAM  STAYL 

002 

C  REVISED  01/28/84 

003 

C  THIS  PROGRAM  COMPUTES  ELEMENT  EXCITATIONS  FOR 

TAYLOR 

004 

C  TYPE  SIDELOBES  USING  SYNTHESIS  EXPRESSIONS  OF 

SHELTON 

005 

DIMENSION  Z(  1 00 ) , AM( 1 00 ) , EN( 1 00 ) , EX(  1 00 ) 

006 

WRITE  (1,100) 

007 

100 

FORMAT  (•  ENTER  N,NBAR,SIDELOBE  LEVEL  FOR 

TAYLOR* 

008 

1  ,  ’  SYNTHESIS’  ) 

009 

READ  (1,200)  N,  NBAR,  SLL 

010 

200 

FORMAT  (2I0.F0.0) 

011 

WRITE  (4,300)  N,  NBAR,  SLL 

012 

300 

FORMAT  (’  TAYLOR  SYNTHESIS  -  SHELTON’  /  ’ 

N=*,I5, 

013 

1  2X , ’ NBAR= ’ , I 5 , 2X , ’ SIDEL0BE= ' ,  F5.2) 

014 

AL2  =  0.30102999566398 

015 

ALE  =  0.43429448190325 

016 

PI  =  3.14159265358979 

017 

M  =  (N-1)/2  +  0.1 

018 

IE  =  1 

019 

IF  (N  ,EQ.  (2*M+1 ) )  IE  =  0 

020 

A  =  (SLL  +  20.0*AL2)  /  (20.0*PI*ALE) 

021 

XN  =  FLOAT  (N) 

022 

XN12  =  FLOAT  (NBAR)  -  0.5 

023 

N1  =  NBAR  -  1 

024 

ALPHA  =  SORT  (A«A  +  XN12  •  XN12) 

025 

DO  1  1=1, N1 

026 

XI 12  =  FLOAT  (I)  -  0.5 

027 

BETA  =  SQRT  (A*A  +  XII 2  *  XII 2) 

028 

Z(I)  =  ( (2.0*PI/XN)/ALPHA)  «  FLOAT  (NBAR) 

»  BETA 

029 

1 

CONTINUE 

030 

DO  2  I=NBAR,M, 1 

031 

Z(I)  =  (  FLOAT  (I)  *  2.0*PI)/XN 

032 

2 

CONTINUE 

033 

EO  =  1.0 

034 

DO  3  1=1, M 

035 

3 

EO  =  EO  «  (1 .0-C0S(Z(I))) 

036 

DO  5  1=1, N1 

037 

AM(I)  =  1.0 

038 

DELTA  =  (2.0*PI  *  FLOAT(I) )/XN 

039 

IF  (IE  .EQ.  1)  AM(I)  =  COS  (DELTA/2.0) 

040 

DO  4  J=1 ,M 

041 

4 

AM(I)  =  AM(I)  »  (COS (DELTA)  -  COS(Z(J))) 

042 

AM(I)  =  AM(I)/EO 

043 

5 

CONTINUE 

044 

DO  6  1=1, M+1 

045 

XI  a  2«I  -  2 

046 

IF  (IE  .EQ.  1)  XI  s  XI  ♦  1 

047 

EN(I)  =  0.0 

048 

DO  7  J=1 ,N1 

049 

XJ  =  FLOAT  (J) 

050 

EN(I)  =  AM( J)  •  COS( (PI*XI*XJ)/XN)  +  EN(I) 

051 

7 

CONTINUE 

052 

EN(I)  =  2.0  •  EN(I)  +  1.0 

053 

6 

CONTINUE 

054 

DO  50  K=1 ,M+1 
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055 

L  =  N+1-K 

056 

EX(K)  =  EN(M+2-K)/EN(M+1) 

057 

50 

EX(L)  =  EX(K) 

058 

WRITE  (4,301) 

059 

WRITE  (4,55)  (I,EX(I) ,  1=1, N) 

060 

301 

FORMAT  (»  ELEM.  NO.*,  3X,  'EXCITATION') 

061 

55 

FORMAT  (5X,I2,5X,F14.7) 

062 

STOP 

063 

END 

J 
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*  1  *■ 


001 

c 

PROGRAM  CONTAYL 

002 

C«<H 

IH#« 

003 

c« 

LINEAR  ARRAY  SYNTHESIS  USING  CONVOLUTION  METHOD. 

004 

c» 

TAYLOR  SIDELOBE  DESIGN 

005 

c» 

ODD  NUMBER  OF  ELEMENTS. 

006 

C«HH 

IIHt 

007 

REAL  PSI(IOO) 

008 

REAL  C(100),AA(100),A1(100),A2(100),A3(100),C0NVO00) 

009 

DATA  A1/100*1 ./,  A3/100*1./,  CONV/100*0./ ,  AA/100*0./ 

010 

Oil 

c» 

N  =  NUMBER  OF  ELEMENTS  IN  THE  ARRAY.  MUST  BE  ODD!! 

012 

c» 

SLL  =  SIDE  LOBE  LEVEL  IN  DBS. 

013 

I#### 

014 

PI  =  3.1415297 

015 

WRITE  (1,100) 

016 

100 

FORMAT  ('  ENTER  DATA:  N , NBAR , SLL  * ) 

017 

READ  (1,200)  N, NBAR, SLL 

018 

200 

FORMAT  (2I0,F0.0) 

019 

M  =  N-1 

020 

NN  =  (N+O/2 

021 

MM  =  (N-D/2 

022 

MD  =  (MM/2)  +  1 

023 

10 

CALL  TAYLX ( MM , SLL , PI , PSI , N , NN , M , NBAR ) 

024 

I##*## 

025 

c* 

CREATE  THREE  ELEMENT  ARRAYS 

026 

Iff###* 

027 

30 

DO  40  I  :  1 ,MM 

028 

A2(I)  =  -2.  •  COS(PSKI)) 

029 

40 

CONTINUE 

030 

c*»< 

>##*•# 

031 

C* 

REPEATED  CONVOLUTION  OF  3-ELEMENT  ARRAYS 

032 

c»*< 

imi* 

033 

CONV(1)  =  A1(1) 

034 

CONV(2)  =  A2( 1 ) 

035 

CONV( 3)  =  A3(1) 

036 

L  *  1 

037 

K  =  5 

038 

LX  =  0 

039 

50 

LL  =  NN-L 

040 

LX  =  LX+1 

041 

60 

L  =  LL 

042 

C(1)  =  A1(L) 

043 

C(2)  =  A2(L) 

044 

C(3)  *  A3(L) 

045 

DO  70  I  =  1 ,3 

046 

DO  70  J  =  I,K 

047 

JJ  =  J-I+1 

048 

AA(  J)  =  AA( J)  +  CONV(JJ)»C(I) 

049 

70 

CONTINUE 

050 

DO  80  I  =  1 ,K 

051 

CONV(I)  =  AA(I) 

052 

AA(I)  a  0.0 

053 

80 

CONTINUE 

054 

K  a  K+2 

C55 
C56 
057 
058 
059 
060 
061 
062 
063 
06  4 
065 
066 
067 
068 
06? 
070 
071 
072 
073 
074 
075 
076 
077 
078 
079 
080 
081 
082 
083 
084 
085 
086 
087 
088 
089 
090 
091 
092 
093 
094 
095 
096 
097 
098 
099 
100 
101 
102 

103 

104 

105 

106 

107 

108 


IF  (L.EQ.MD)  GO  TO  90 
LL  =  NN+1-L 
LSUM  =  L+LX 

IF  (LSUM.EQ.NN)  GO  TO  60 
GO  TO  50 
90  CONTINUE 

C  WRITE  (4, KOI) 

C  401  FORMAT  (  '  PSI  ZEROS'/) 

C  DO  110  I  =  1 ,M 

C  WRITE  (4,500)  I,PSI(I) 

C  110  CONTINUE 

500  FORMAT  (10X,I2,10X,F10.6) 

WRITE  (4,600) 

600  FORMAT  (  •  CURRENTS'/) 

DO  120  I  =  1  ,N 

WRITE  (4,700)  I , CONV ( I ) 

120  CONTINUE 

700  FORMAT  ( 1 0X, 12, 1 0X,F10. 6) 

STOP 

END 

FUNCTION  COSH(R) 

Y  =  EXP(R) 

COSH  =  (Y  +  ( 1 ,0/Y) )/2. 

RETURN 

END 

C«  INVERSE  HYPERBOLIC  COSINE  FUNCTION 
C»**»**»* 

FUNCTION  RCOSH(R) 

RCOSH  =  ALOG(R  +  SQRT(R*R  -  1.0)) 

RETURN 

END 

C*  COMPUTATION  OF  TAYLOR  ZEROS 

C**”»«* 

SUBROUTINE  TAYLX(MM , SLL , PI , ZEROS , N , NN , M , NBAR ) 

REAL  ZER0(50),  ZEROS(IGO),  MEMA , MEMB 
NBAR1  =  NBAR-1 
A  =  (SLL  +  6.0202)/27.2875 

C»  COMPUTE  ZEROS  FROM  1  TO  NBAR 

DO  10  I  =  1 ,NBAR1 
RI  =  I 

MEMA  =  (A»A)  +  ( (RI-.5)**2.) 

MEMB  =  (A«A)  +  ( (NBAR-.5)**2. ) 

ZERO( I)=( ( (2.*PI)*NBAR)/N)*( ( SQRT( MEMA ) ) / ( SQRT( MEMB ) ) ) 
10  CONTINUE 

C»  COMPUTE  ZEROS  FROM  NBAR  TO  M 


DO  20  I  s  NBAR, MM 
RI  =  I 

ZERO(I)  =  (2«PI*RI)/N 


20 


1 


109 

20 

CONTINUE 

110 

C 

WRITE  (4,100) 

111 

C 

100 

FORMAT  (  *  ZEROS') 

112 

C 

DO  30  I  =  1 ,MM 

113 

C 

WRITE  (4,200)  I,ZERO(I) 

114 

C 

30 

CONTINUE 

115 

C 

200 

FORMAT  ( 10X,I2, 10X,F10.6) 

116 

DO  40  J  =  1 ,NN 

117 

ZER0S(NN-1+J)  =  ZERO(J) 

118 

ZEROS(NN-J)  =  ZERO(J) 

119 

40 

CONTINUE 

120 

RETURN 

121 

END 
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